Attractive vortex interaction and the intermediate-mixed state of superconductors 
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The magnetic vortices in superconductors usually repel each other. Several cases are discussed 
when the vortex interaction has an attractive tail and thus a minimum, leading to vortex clusters 
and chains. Decoration pictures then typically look like in the intermediate state of type-I super- 
conductors, showing lamellae or islands of Meissner state or surrounded by Meissner state, but with 
the normal regions filled with Abrikosov vortices that are typical for type-II superconductors in the 
mixed state. Such intermediate-mixed state was observed and investigated in detail in pure Nb, TaN 
and other materials 40 years ago; last year it was possibly also observed in MgB2 , where it was called 
"a totally new state" and ascribed to the existence of two superconducting electron bands, one of 
type-I and one of type-II. The complicated electronic structure of MgB2 and its consequences for su- 
perconductivity and vortices are discussed. It is shown that for the real superconductor MgB2 which 
possesses a single transition temperature, the assumption of two independent order parameters with 
separate penetration depths and separate coherence lengths is unphysical. 

PACS numbers: 



1. VORTEX LATTICE FROM GL THEORY 



The most successful phenomenological theory of su- 
perconductors was conceived in 1950 by Ginzburg and 
Landau (GL) [l|. When written in reduced units (length 
unit A, magnetic field unit \2H C , energy density unit 
HoHc, where A is the magnetic penetration depth and 
H c the thermodynamic critical field) the GL theory con- 
tains only one parameter, the GL parameter re = A/£ (£ 
is the superconducting coherence length). The kind of 
solutions of GL theory is quite different when re < l/v2 
or re > l/v2. Physically, this is due to the fact that 
the energy of the wall between a normal and a super- 
conducting domain is positive when re < 0.71 (type-I) 
and negative when re > 0.71 (type-II). This implies that 
a type-II superconductor is unstable against the sponta- 
neous formation of many small superconducting and nor- 
mal domains. A thesis student of Lev Landau in Moscow, 
Alexei Abrikosov in 1953 discovered a periodic solution 
of the GL theory which describes the occurrence of a reg- 
ular lattice of vortex lines (published in 1957 |2|). Each 
Abrikosov vortex (or fluxon, flux line) carries one quan- 
tum of magnetic flux $ = h/2e = 2.07 • 10~ 15 Tin 2 , 
and the supercurrent circulates around a singular line on 
which the complex GL function ip(x,y,z) (or supercon- 
ducting order parameter l^l 2 ) is zero. For finding this 
vortex lattice solution, in which Landau initially did not 
believe |2[ , Abrikosov obtained the Nobel prize in physics 
50 years later in 2003. 

Within the original GL theory in bulk superconduc- 
tors, with increasing applied magnetic field B a = /ioH a 
the superconductor is first in the Meissner state which 
has induction B — 0. When B a reaches the lower critical 
field B c \ Abrikosov vortices start to penetrate. With fur- 
ther increasing B a the vortex lines form a more and more 



dense triangular lattice with induction B = (2/^/3)^> /a 2 
(a = vortex spacing). When B a and B reach the upper 
critical field B C 2 = fJ.oH c2 = ^o/{2'n'£, 2 ), the order param- 
eter vanishes and the superconductor becomes normal 
conducting. Since B(B a ) increases monotonically (see 
Fig. 1 (top)), one may say that the pressure of the vortex 
lattice is positive and the vortices repel each other, held 
together by the applied magnetic field B a , which acts as 
an external pressure. From thermodynamics follows that 
in equilibrium (in absence of pinning and surface bar- 
riers) one has B a = HoH a = fi Q dF/dB where F(B, re) 
is the GL free energy density. The numerical GL solu- 
tion B a (B) of the vortex lattice Q for 1/^/2 < re < 10 
may be fitted well by the expression (with b = B/B c2l 
h = B a /B c2l m = b-h = M/B c2 , h cl = B cl /B c2 ): 
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with (3a = 1.15960 (1.18034) for the triangular (square) 
vortex lattice. At re = l/v2 the B(B a ) is a vertical 
straight line at B a = B c \ — B c2 like in type-I supercon- 
ductors. This means the vortices do not interact with 
each other and all vortex arrangements have the same 
free energy, when re = l/v2- 

The vanishing vortex interaction for re = l/v2 can 
also be seen from Kramer's [4| asymptotic interaction 
between two vortices at large distance r 3> A, V(r) = 
di(K)Ko(r/X) — d 2 (K)K (V2r/t;) where the constants are 
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FIG. 1: Top: Ideal reversible magnetization curves of a type- 
II superconductor with no vortex pinning and no demagne- 
tizing effects (long cylinder or slab in a parallel magnetic field 
B a ) from Ginzburg-Landau Theory. Plotted is the induction 
B versus B a for several GL parameters k = l/\/2 ... 3. Bot- 
tom: The same magnetization curve B(B a ) from the micro- 
scopic BCS-Gor'kov theory for GL parameter k « 0.85 and 
temperature T/T c ~ 0.7, schematic from [281 ] . The dashed 
part for B < Bo is unphysical, the real B(B a ) contains the 
depicted vertical line from B — to B = Bo obtained from a 
Maxwell construction that divides the area under the dashed 
curve into equal halves. 




FIG. 2: Intermediate state of a type-I superconductor. The 
normal domains are dark. Coexistence of triangular lattice of 
flux tubes and of laminar domains. Tantalum disk of tickness 
d = 33 fim, diameter D = 5 mm, T = 1.2 K , B a = 34 mT. 
(Courtesy U. Essmann). Optical microscope. 
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d\ = d,2 at k — 1/V2. A similar vortex potential [5( 
was derived such that it reproduces the nonlocal elastic 
properties of the triangular vortex lattice [6( in the entire 
ranges of B and k and for all wavelengths of the elastic 
strain. For parallel vortices this approximate interaction 
potential reads 



V(r) = 



& 



2ir\' 2 fi 



[K a (r/\i)-K (r/O\ 



(2) 



wit h the len gths A' = X/(\^\ 2 ) 1/2 « A/v 7 ! - !, £,' = 
^/•v/2(l — b) and Kq(x) a modified Bessel function with 
the limits — lnx (x <C 1) and (ir/2x) 1 / 2 exp(— x) (x > 1). 
Physically, the first (repulsive) term in V(r) Eq. (2) de- 
scribes the magnetic interaction and the second (attrac- 
tive) term originates from the overlap of the vortex cores 
that reduces the condensation energy ex 1 — |V>| 2 - Note 
that this V(r) has a finite maximum V(0), approximately 



FIG. 3: Laminar domains in the intermediate state of pure 
superconductors at T — 1.2 K. a) Type-I superconductor, 
Tantalum disc, d = 33 /im, D — 4 mm, B a = 68 mT, angle 
between magnetic field and disc 15°. Optical microscope, b) 
Type-II superconductor, Niobium disc, d = 40 /im, D = 4 
mm, B a — 74 mT. Optical microscope, c) As b) but electron 
microscope. From [11]. 



equal to half the self energy of a vortex with flux 2$ 0j 
since the two terms diverging as ±lnr compensate. 

For k > l/\/2 the potential (2) is repulsive. At k = 
l/-\/2 it vanishes at all distances, which means that all 
vortex arrangements have the same energy. This is an 
exact result of the GL theory [7|. The particular case 
« = l/v2 and its surroundings (k « 0.71, T sw T c ) were 
investigated in detail in [8|, |9( . 
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FIG. 4: Flux-line lattice observed at the surface of type-II 
superconductors in an electron microscope after decoration 
with Fe microcrystallites ("magnetic smoke"). Top: High- 
purity Nb disks 1 mm thick, 4 mm diameter, of different crys- 
tallographic orientations [110] and [Oil], at T = 1.2 K and 
B a — 800 Gauss (B c i = 1400 Gauss). Due to demagnetiz- 
ing effects and the small k « 0.70, round islands of Meissner 
phase are surrounded by a regular vortex lattice ("interme- 
diate mixed state"). Bottom: High-purity Nb foil 0.16 mm 
thick at T = 1.2 K and B a = 173 Gauss. Round islands of 
vortex lattice embedded in a Meissner phase. (Courtesy U. 
Essmann) 



2. TYPE-I SUPERCONDUCTORS 

Superconductors with k < l/\/2 are called as type 
I. For B a < B c = noH c (thermodynamic critical field) 
these should exhibit ideal Meissner effect with B = 
(except for a surface layer of thickness A) and B = B a 
for B a > B c . But in samples of finite size the demagneti- 
zation factor N > allows the observation of an "inter- 
mediate state" containing superconducting and normal 
domains with induction B — and B = B C1 respectively, 
in the range (1 - N)B C < B a < B c [HJ. Figure 2 shows 
a decoration picture of the type-I superconductor Tanta- 
lum, a disk with thickness d — 33 //m, diameter D = 5 




FIG. 5: Decoration of a square disk 5x5x1 mm of high- 
purity polycrystalline Nb at T — 1.2 K and B a — 1100 Gauss. 
As B a is increased, magnetic flux penetrates from the edges 
in form of fingers which are composed of vortex lattice shown 
enlarged in the right picture. The rectangular cross section 
of the disk causes an edge barrier (section 8.5.). As soon as 
this edge barrier is overcome, single flux lines or droplets of 
vortex lattice (lower right) pull apart from these fingers and 
jump to the center, filling the disk with flux from the center. 
(Courtesy U. Essmann) 



1.2 K and B a = 34 mT, taken in 1969, from 
The picture shows the coexistence of 
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triangular lattice of flux tubes and of laminar domains. 
As recently shown by Prozorov 12|, [13| , the equilibrium 
topology of the intermediate state are flux tubes (in ab- 
sence of bulk pinning and geometric barrier). Interest- 
ingly, pin-free disks and strips of constant thickness show 
hysteretic magnetization curves due to edge barriers for 
penetrating flux tubes and exiting lamellae. This topo- 
logical hysteresis of type-I superconductors was observed 
and calculated analytically already in 1974 [14[, where it 
was shown that ellipsoidal samples exhibit no magnetic 
hysteresis. Magnetic hysteresis loops in pin-free type-I 
superconductor disk and strips look very similar as the 
irreversible magnetization loops in pin-free type-II super- 
conductors with edge barrier [15], in particular, they go 
through the origin (M = B — B a = at B a = 0, no 
flux is trapped) and become reversible at large B a (at 
B a > B cl /2 for type-II). 

The structure of the flux tubes or lamellae in type-I 
superconductors and their widening (mushrooming) and 
branching as they approach the surface were calculated, 
e.g., in [lg, |l7[. In thin films of type-I superconductors a 
vortex lattice occurs though k < 0.71 [18|. In particular, 
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FIG. 6: Ideal magnetization curves of long cylinders or slabs 
in a parallel magnetic field B a (demagnetization factor N — 0, 
left) and of a sphere (N = 1/3, right). Top row: Type-I 
superconductor with Meissner state for B a < (1 — N)B C , in- 
termediate state for (1 — N)B C < B a < S c , and normal state 
for B a > B c . The dashed line indicates the enhancement 
of the penetration field when the wall energy of the normal 
and superconducting domains in the intermediate state is ac- 
counted for. Middle row: Type-II/1 superconductors with 
Ginzburg-Landau parameter re close to l/s/2, for N = ex- 
hibit a jump of height Bo in their induction and magnetiza- 
tion at B a = -Bci due to an attractive interaction between flux 
lines. For N > 0, this jump is stretched over a finite range of 
B a , allowing one to observe an intermediate mixed state with 
domains of Meissner phase and vortex-lattice phase. Bottom 
row: Type-II/2 superconductors with re 3> l/v2- Meissner 
phase for B a < (1 — N)B c i, mixed state (vortex lattice) for 
(1 — N)B C \ < B a < B C 2, and normal state for B a > B C 2- 



for k = 0.5 the shear modulus cqq of this vortex lattice 
is positive at all inductions < B < B C 2 when the film 
thickness is d < £, see Figs. 9 and 13 in [l9j, and the 
magnetization curve —m{h) = —m(b) (Fig. 4 in [19j) has 
positive curvature in this case. These GL results show 
that in thin (and infinitely large) films in perpendicular 
field B a — B the vortices may repel each other even when 
K < 0.71 since the long-range repulsive interaction via 
the magnetic stray field dominates [20J. For the vortex 



interaction in thin films of finite size see 21 1 



3. ATTRACTIVE VORTEX INTERACTION 

The first successful decoration experiments evaporated 
an iron wire in a 1 Torr He atmosphere that produced 
"magnetic smoke" which settles on the surface of the su- 
perconductor and marks the ends of the flux tubes or 
flux lamellae (Fig. 2,3a) or of the vortices (Figs. 3c, 4, 5) 



22j. In the type-II superconductors Vanadium and Nio- 
bium (see Fig. 3 of llj) and Pbln alloys they showed 
more or less regular triangular vortex lattice with struc- 
tural defects, e.g., dislocations and vacancies. These pic- 
tures were consistent with the prediction from GL theory 
that the vortices repel each other. However, later deco- 
ration of very pure Nb showed a domain structure simi- 
lar to that observed in type-I superconductors but with 
the normal domains replaced by domains of vortex lat- 
tice. The Figs. 3a and 3b compare the laminar domains 
in the type-I superconductor Tantalum and in the type- 
II superconductor Nb. In an optical microscope these 
domains look very similar, but the electron microscope 
(Fig. 3c) reveals that in Nb the domains which carry the 
magnetic flux consist of vortex lattice with constant lat- 
tice spacing clq. 

In Fig. 4 (top) one sees islands of Meissner state em- 
bedded in a vortex lattice, while Fig. 4 (bottom) shows 
Meissner state with islands of vortex lattice. In Fig. 5 
tongues of vortex lattice penetrate from the specimen 
edge into the Meissner state of high-purity polycrystallinc 
Nb. 

This observation of an upper limit clq of the vortex 
spacing in the intermediate state of low-re type-II super- 
conductors, and its explanation by an attractive inter- 
action, suggests that the magnetization curves in the 
ideal (pin-free) case should exhibit a jump of height 
B = (2/v / 3) ( f > o/ a o at tne ne ld of first vortex penetra- 
tion. This means, as soon as it is energetically favorable 
for vortices to penetrate, they immediately jump to their 
equilibrium distance ao, which for N = means a uni- 
form induction Bq, see Fig. 1 (bottom). Such jumps of 
B(B a ) from to Bq were indeed observed in clean Nb 
and in lead alloys with increasing content of Thallium 
(which changes k), and excellent agreement of Bq with 
the directly measured a was observed 23H25|. Jumps 
Bq, and thus vortex attraction, were also observed in the 
TaN system with varying nitrogen content and various 
k = 0.35 to k = 1.05 Q. See also the reviews QHIlll. 

The jump B in the ideal magnetization curve B(B a ) 
was confirmed by computations [28( based on the mi- 
croscopic BCS-Gor'kov theory that reduces to the GL 
theory in the limit T — >• T c . Earlier computations [29j 
confirming the vortex attraction were based on the ex- 
tended GL theory of Neumann and Tewordt [30| that 
keeps all correction terms to lowest order in 1 — T/T c . 
The BCS-Gor'kov computations [28j for clean supercon- 
ductors yield an S-shaped curve B(B a ) with an unstable 
part (dashed curve in Fig. 1 (bottom)) that has to be 
replaced by a vertical line of height Bq that cuts the area 
under the dashed line in two equal parts (Maxwell con- 
struction). Useful interpolation formulae for this S-shape 



are given in [31 1 



The phase diagram of low-re clean superconductors in 
the K—t plane (t = T/T c ) was measured by Auer and Ull- 
maier [26( and computed from BCS-Gor'kov-Eilenberger 



theory by Klein [32j; this is discussed in [33|. Near 
K = 0.71 with increasing t, this phase diagram shows the 
transition from type-I to usual type-II (termed type-II/2) 
via a region where the vortices attract each other, termed 
type-II/1 superconductors. Figure 6 shows the magneti- 
zation curves —M(B a ), M = B—B a , for type-I, 11/ 1, and 
II/2 superconductors shaped as cylinders (N — 0, left 
column) and spheres (N = 1/3, right column), schematic 
from [27j . Recently by neutron-scattering a kind of phase 
diagram of vortex lattices and normal and intermediate- 
mixed states has been obtained for pure Nb cylinders of 
14 mm length and 4 mm diameter 3_4|. So, vortex at- 
traction in pure Niobium is still an active field. 

The above attractive interaction originates from BCS- 
corrections to the GL theory (which, strictly spoken, is 
valid only at T — T c ) and at < T < T c occurs for 
0.71 < k < 1.5 [11]. The London theory that follows 
from GL theory when k 2 ^> 1 and b <C 1, has purely 
repulsive vortex interaction, namely, the first term in 
Eq. (2). However, many corrections to the simple GL 
or London theories are expected to modify the monoton- 
ically decreasing interaction potential at large distances, 
V(r) oc exp(— r/A), such that A becomes complex. This, 
in principle, causes an oscillating potential, whose first 
minimum may occur at large distances where the ampli- 
tude of the potential is small. 

Such corrections may have various origins, e.g., the 
spatially varying order parameter may lead to non-local 
electrodynamics [35j. Linear nonlocal electrodynamics 
relates the supercurrent density j to the vector poten- 
tial A (in Fourier space and London gauge) by a kernel 
Q(k), j(k) = -Q(fc)A(k). The simple "local" London 
theory has Q(k) = Q(0) = A~ 2 . An elegant non-local ex- 
tension of London theory uses the Pippard kernel Qp(k) 
36], which introduces a Pippard length £p. This Pip- 
pard theory leads to a vortex field H(r) and vortex po- 
tential V(r) = H(r)<&o that oscillates at large distances 



( "field reversal" |37j,|38|). The corresponding BCS kernel 
QBCs(k) is an infinite sum of such Pippard kernels and 
its vortex field H{r) and potential V(r) exhibit similar 
oscillations [37j , see also 13811 . We note that some of the 
early theories quoted in [261 ] were flawed. 

In anisotropic superconductors vortex attraction may 
even occur in the linear London theory, which now has 
different penetration depths A a , At,, and A c for supercur- 
rents flowing along the crystalline a, b, and c axes [27[. A 
field reversal and vortex attraction may then occur when 
the vortices are at an oblique angle [39H41]. 

In layered superconductors an anisotropic pair po- 
tential between the pancake vortices is investigated by 
Feigelman et al. 42j and Clem [43J. This potential has 
an intra-layer logarithmic repulsive part and an inter- 
layer attractive part, both arising from the magnetic in- 
teraction. This potential is used by Jackson and Das 
[44 |45| for the study of freezing of vortex fluids by a 
density functional method, yielding a nearly first order 



transition from a fluid to a triangular Abrikosov vor- 
tex lattice. Brandt et al. [46J predicted long-ranged 
fluctuation-induced attraction of vortices to the surface 
in layered superconductors by a type of Casimir force. 
This work inspired Blatter and Geshkenbein J47J to study 
long-range van der Waals type inter-vortex attraction 
in the anisotropic layered materials, which arises from 
the thermal fluctuations of vortices. At low magnetic 
fields, the inter-vortex separation being very large, even 
weak attraction influences the phase diagram near the 
H c i boundary. While the above forces are of thermal 
origin, van der Waals attraction mediated by impurities 
is studied by Mukherji and Nattermann [48j. 

It is conceivable that also a two-component GL theory 
may lead to an oscillating vortex- vortex interaction with 
attractive tail and more or less pronounced minimum if 
its input parameters are chosen appropriately. Such a po- 
tential was calculated in [49J, |50{ ■ In the same work some 
of the decoration [49] and scanning SQUID microscopy 

50] images taken during field cooling of MgB 2 platelets 
suggest that at the (unknown) temperature where the de- 
picted vortex positions during field cooling were frozen, 
the vortex interaction had a minimum that determined 
the dominating vortex distance in the observed vortex 
chains and clusters. 

While Refs. 49j, |50j do not mention any of the above 
listed other cases of vortex attraction, the supplement 

51| refers to previously observed and well understood 
vortex attraction in pure Nb. Note that possible obser- 
vation of vortex attraction in MgB2 does not mean that 
the input parameters and interaction potential of simu- 
lations that lead to similar vortex arrangements are true 
or unique. In the following the case of the 2-band super- 
conductor MgB2 will be discussed in some detail. 



4. BEYOND THE STANDARD BCS MODEL 

The BCS theory is the standard model for understand- 
ing various basic properties of conventional metallic su- 
perconductors. In this simple theory the Fermi surface is 
isotropic. The attractive electron-phonon interaction is a 
constant in a small shell around the FS. In 1958 Gorkov 
showed the equivalence of the BCS theory with the one- 
order parameter isotropic Ginzburg-Landau theory. The 
latter is a powerful phenomenological theory suitable for 
practical applications. On the other hand a more rigor- 
ous approach is the Eliashberg theory suitable for strong 
coupling of electron-phonon interaction. 

In the past three decades a number of novel super- 
conductors have been discovered with many exotic prop- 
erties beginning with the high-T c cuprates. Others are 
doped fullerides, ruthenides, nickel borocarbides, magne- 
sium diborides and pnictides etc. Basic understanding of 
superconductors requires the correct information about 
its electronic and vibrational structures and electron- 



phonon interaction. Fortunately an ab-initio theory like 
the density functional theory helps us to meet this type of 
requirement. By using a self-consistent elctron potential 
arising out of ionic, Hartree, exchange and correlations, 
one computes the quasi-particle dispersion throughout 
the Brilliouin zone and obtains the Fermi surface. 

In a superconductor the formation of Cooper pairs de- 
pends heavily on the geometry of the Fermi surfaces. 
Then comes the electron-phonon interaction that medi- 
ates the pairing. In view of the complicated band struc- 
tures and their Fermi surfaces, the novel superconductors 
are treated beyond the standard BCS model. We are 
reminded that single-sheeted Fermi surfaces only come 
from the simple metals like alkalies considered in the jel- 
lium model and those are not superconductors under am- 
bient conditions. Therefore, if the multi-sheeted Fermi 
surfaces are connected in the entire Brilliouin zone, then 
an anisotropic GL theory is useful. Otherwise for uncon- 
nected Fermi surfaces one can use either (i) a multi-band 
BCS theory 52[ or (ii) a generalization of the GL theory 
invoking many (or many components of) order param- 
eters [53j- This two-component order- parameter theory 
was used to explain successfully two specific-heat transi- 
tions in the high T c superconductor YBCO (see Choy et 
al. in (H). 

As said before, in reality the bands are formed due 
to ionic and electronic self-consistent potential. In prin- 
ciple we are dealing with a charged (Coulomb) system, 
hence the bands are not entirely independent, and further 
due to the electron-phonon interaction (whether small 
or large) the bands are coupled. For this reason multi- 
(inter- and intra-)band interaction terms explicitly ap- 
pear in the microscopic Hamiltonian |52j. In the case 
of multi-order parameter GL theory, one has the super- 
conducting free energy composed of individual parts of 
the GL free energy and the additional free energy due 
to coupling between the order parameters as mentioned 
above. This choice of coupling is non-unique. If ipx and 
t/>2 are two order parameters the coupling or interaction 
term can be one or more terms from the following list: 

iMi\ 2 m\ 72(n*v+nv> 2 + hc), 73(^2 + ^2") 

with n = V + 27riA/$ - The first term of this list is 
a systematic expansion of GL, the second term means 
coupling through the gradient of order parameter and 
vector potential A, and the third term is the internal 
Josephson coupling term. This last term is supposed to 
provide a "minimal coupling" [54|-l56| . Instead of the phe- 
nomenological GL theory, we start with the microscopic 
Bogolyubov pairing Hamiltonian and use the mean-field 
approach; we derive the GL functional with a coupling 
term that will be Josephson-like. An important point 
here is that the coefficients a, of \ipi\ , and the coefficient 
73 of the Josephson term, do depend both on the intra- 
band and inter-band pairing. Therefore, as a matter of 
principle the order parameters can never be independent. 



There is a huge amount of literature on both microscopic 
and GL theory for the multi-band cases. We do not in- 
tend to review that here, see for example[55j. We shall 
return to the interacting bands in the next Sections. 

In the earlier sections we have discussed the physics of 
vortex matter with particular emphasis on the value of k 
and the type of superconductivity (type-I or type-II) . We 
also discussed the behavior of pure Nb at temperatures 
below the GL regime. In the conventional single gap 
superconductor within GL theory there are a well defined 
coherence length £ and penetration depth A. In terms of 
them one defines k = A/£. As stated earlier its unique 
value, if less or bigger than 0.71, defines precisely type- 
I or type-II superconductors, in which the surface free 
energy between normal and Meissner domains is positive 
or negative, respectively. This definition tells about the 
equilibrium states of a superconductor. Here with this 
brief background we shall discuss the interesting physics 
of the novel superconductor MgB 2 and the vortex matter 
therein. Needless to say that it has become a news item 
in the past one year with the claim that a new form of 
matter is discovered 49|, [50j . Detailed discussions follow 
in Sec. 6. 



5. VORTICES IN MGB 2 

In 2001 by a surprising discovery Nagamatsu et al. [57j 
reported T~ of an intermetallic MgB2 around 39 K. Soon 
after (see [58[ database for the first three years) it was 
realized that MgB2 offers a host of novel properties that 
has been explored during the past 9 years. Magnesium 
diboride has a crystallographic structure consisting of 
honeycomb boron layers separated by magnesium layers, 
which have also honeycomb structure. Magnesium atoms 
are ionized in this structure donating their s-electrons to 
the conduction band. There are two 3-dimensional 7r 
bands arising from the boron p z orbitals. One of them 
is electron-like, whereas the other is hole-like. From the 
Px,y orbitals two a bands occur. They are 2-dimensional 
confined to the boron planes. Bonds within the boron 
layers are strongly covalent, whereas between the layers 
they are metallic. From this a complicated Fermi surfaces 
geometry emerges. Two cylinders around the T-A-r lines 
are two a bands, and two webbed tunnels are attributed 
to 7r bands, see Fig. 7 (from Mazin and Antropov [59(). 
This overall picture is consistent with all the electronic 
structure calculations 59]. 

It is further found [60( that the optical bond-stretching 
Ei g phonons couple strongly to the holes at the top 
of a bands, whereas the 3-dimensional 7r electrons 
are weakly coupled to the phonons. For this reason 
of different coupling strengths of a and 7r bands, a 
two-gap superconductor merges. The E2 g phonons 
strongly coupled with a band, produce high T c in 
MgB2. A further coupling with the 7r bands enhances 




FIG. 7: Band structure and Fermi surface of MgB2 from 
Mazin and Antropov [K|. Green and blue cylinders (hole- 
like) are from the a bands, and the blue (hole-like) and red 
(electron-like) tubular networks are from the 7r bands. 



T c . Table 1 presents values of intra- and inter-band 
(electron-phonon) coupling-constant matrix-elements 
calculated by various groups. The numbers show that 
the electron-phonon coupling A (a dimensionless param- 
eter [6l| not to be confused with the penetration depth 
A) for the a band is very strong and for the 7r band 
it is weak. The inter-band coupling is still weaker but 
appreciable. These numbers are useful ingredients to be 
incorporated in the two-order parameter calculations. 
In a first principles calculation using Eliashberg theory, 
Geerk et al. |60| emphasize the importance of inter-band 
pairing and low energy optical modes arising out of 
boron atoms vibrating perpendicular to the boron planes. 



A„ 



At, 



Xn 



A71 



1.017 0.448 0.213 0.155 Golubov et al. [60] 
0.96 0.29 0.23 0.17 Liu et al. [60] 
0.78 0.21 0.15 0.11 Choi et al. [59] 

Table 1. Matrix elements of electron-phonon interactions 



In the past nine years a variety of experiments have 
shown that MgB 2 is a two band superconductor. The 
main experiments are specific heat (Choi et al. i n 1591 l62j . 
point contact and tunnelling spectroscopy [63J, l64j , Ra- 
man spectroscopy [651 ] . penetration depth measurements 
[66j, angle-resolved photo-emission spectroscopy [671 ] . 
small-angle neutron-scattering J68| . and muon-spin re- 
laxation studies [69| etc. The value of the gap on the 
a bands A CT ranges from 6.4 to 7.2 meV, while on the 7r 
bands A„- varies from 1.2 to 3.7 meV. Almost all param- 
eters such as coherence lengths (£ CT a nd £ „■) [661 ] . upper 
critical fields (H c2 (a) and H r oM) [70, [71], and penetra- 
tion depths (A CT and A ff ) [66j are anisotropic. Besides 
this the anisotropic factor is also strongly temperature 
dependent. These experiments convincingly suggest that 
MgB2 is a two band superconductor, where phonons me- 
diate pairing in the respective bands. Further, it is in- 
vestigated that inter-band coupling helps to increase T c . 
In the absence of inter-band coupling, it may give two 
isolated superconductors with two lower T c . 



6. A NEW PHASE OF MGB 2 ! 

In a recent paper Babaev and Speight [72} studied 
phenomenologically a two-band superconductor with two 
order parameters, starting with a free energy that in- 
cludes a Josephson-like interaction 73 between the two 
bands. Choosing 73 — » in the decoupled band limit, 
they considered Ki (i — 1,2) in two different regimes 
to produce type-I («i < 0.71) and type-II (n\ < 0.71) 
materials. Their prediction leads to, what they call a 
"semi-Meissner state", particularly in the case of very 
low magnetic fields. Babaev and Speight claim that this 
state is thermodynamically stable. Instead of homoge- 
neous distribution, the vortices form aperiodic clusters or 
vortexless Meissner domains, arising out of short range 
repulsion and long range attraction. Below we shall argue 
that in the particular case of MgB2, this is an unrealistic 
construction. 

In what is claimed as a remarkable manifestation of 
a new type of superconductivity in MgB2, Moshchalkov 
and co-workers [49j, [50| reported last year a novel "type- 
1.5 superconductor", contrary to type-I and type-II su- 
perconductors. It appears these authors are influenced by 
the idea of "semi-Meissner state" of Babaev and Speight. 
Apparently both names carry similar meaning. Let us 
analyse carefully the results of Moshchalkov et al. 

Similar to Babaev and Speight, Moshchalkov et al. use 
a two-band GL functional with an inter-band coupling 
term. They constructed the GL parameters «;, = Aj/£j 
(i stands for a and 7r bands of MgB2) from the mea- 
sured values of energy gaps, Fermi velocities and plasma 
frequencies of respective bands. The values are so cho- 
sen that «:<, > 0.71 (=3.68) and k t < 0.71 (= 0.66). 



Thus, MgB2 has properties of both type-I and type-II 
simultaneously, unlike in the case of Nb where k is close 
to 0.71 as discussed above. This construction is in ab- 
sence of the coupling term, so the bands are fully in- 
dependent/noninteracting. We have stated earlier that 
energy gaps A CT range from 6.4 to 7.2 meV, while the A T 
varies from 1.2 to 3.7 meV. Moshchalkov et al. chose the 
value of An- = 7.1 meV and A T = 2.2 meV from some 
experiments and computed £* from the one-band BCS 
relation £j = hv p / '(2tt 2 Aj) . Then they used calculated 
plasma frequencies of respective bands, from which they 
obtained the London penetration depths A CT (0) = 48.8 
nm and A„-(0) = 33.6 nm and finally «, = Xi/£,i- 
Here we note a few points. 

A. (1) The numerical estimates for £, are obtained by 
using the one-band BCS formula. (2) One band calcu- 
lated plasma frequency is used to compute Aj. One does 
not know how these are relevant to interacting two-band 
MgB 2 . 

B. For critical fields Moshchalkov et al. gave their mea- 
sured values for H cl (0) = 241 mT, H c2 {0) = 5.1 T for the 
a band. From these data they also obtained A CT and £ CT 
from the one band formula. The value of H c i(0) = 241 
mT measured for their sample is much larger than the 
values reported by others, which are in the range of 25 
- 48 mT, see Fig. 25 of Ref. |zy. Another point of con- 
cern is the reported value of thermodynamic critical field 
H c (0) = 230 mT (Z2J], whereas Moshchalkov et al. have 
Hci(0) = 241 mT higher than H c (0). 

C. The A that Moshchalkov et al. estimated are smaller 
by a large factor as can be verified from the Table 1 of 
Golubov et al. [3]. This Table lists the values of A(0) 
from a variety of experiments. Smaller A can make the k 
smaller. This possibility has made h* < 0.71. Another 
important point is that using two A for a single flux line 
is questionable. One has to solve a single equation for the 
magnetic field (in two-order parameter GL/BCS theory) 
to obtain a single A, even in the decoupled band limit. 
One has A oc (J2i \^i\ 2 l rn i)~ 1 ^ 2 i containing a sum of two 
superfluid densities. 

An often heard argument is that whether there is one 
K or two k, there will be only one surface free energy 
between Meissner state and normal state, which would 
be either positive, zero or negative. Thus, supercon- 
ductors are either of type-I or type-II. However, this ar- 
gument tacitly assumes that boundaries occur only be- 
tween Meissner and normal states. But as figures of the 
intermediate-mixed state show (see Figs. 3 to 5), there 
the boundary occurs between Meissner and vortex states. 
In these experiments the wall energy between Meissner 
state and normal state is negative (leading to a vortex 
state) but the wall energy between Meissner state and 
the vortex state is positive, causing the observed domain 
structure. So, there are indeed experimental situations 
that exhibit features of both type-I and type-II behavior. 
As stated in Sec. 3, one condition for this intermediate- 



mixed state to occur is that the vortex interaction has 
an attractive tail, or that the theoretical magnetization 
curve is S-shaped, see Fig. 1 (bottom). But in addition, 
the demagnetization factor of the decorated specimen has 
to be N > 0, else only Meissner state or a uniform vortex 
state will occur; the detailed appearance further depends 
on pinning, on the edge barrier, and on the magnetic 
history. This intermediate mixed state has been known 
for long time from experiments [111, |22j-|26j, |3J, l75| and 
computations [H, [H, EH!!] for pure Nb, TaN, Pbln, 



and PbTl. In this respect, the decoration and scanning 
SQUID results of 49j, [50j are not so new. 

As discussed above, several origins of vortex attrac- 
tion are conceivable, depending on the material, its pu- 
rity, magnetic history, and temperature. Computer sim- 
ulations using such a vortex interaction (short-range re- 
pulsive, long-range attractive) yield vortex arrangements 
that are similar to the observed ones [49j, [5CJ . The details 
of the resulting clusters or chains also depend on the type 
and strength of vortex pinning that is always present and 
should be included in such simulations, in particular at 
very low induction B. 



7. SUMMARY AND CONCLUSIONS 

Type-I superconductors with demagnetization factor 
N > in the field range (1 - N)B C < B a < B c are in the 
intermediate state showing normal conducting lamellae 
or tubes inside which B = B c , embedded in Meissner 
state with B = 0. Type-II superconductors in the field 
range (1 — N)B C \ < B a < B c2 contain vortices which 
repel each other (mixed state) . Under several conditions 
the vortex interaction may have an attractive tail, leading 
to observation of vortex clusters or lamellae similar as in 
type-I superconductors, but filled with vortices, and to 
vortex chains. This "intermediate-mixed state" has been 
known since forty years in pure Nb, TaN [26], Pbln, and 
PbTl [75j at temperatures sufficiently low such that GL 
theory does not apply. The recent possible observation of 
vortex attraction in the two-band superconductor MgB2, 
if one agrees with this interpretation of the decoration 
[49J | and scanning SQUID [50( images, may be due to the 
complexity of the required microscopic theory. But we 
question if these images hint at "a totally new state" that 
"behaves in an extremely unusual way" [49]. It should 
also be considered that at the very low inductions of these 
images even very weak pinning will influence the vortex 
positions and may cause vortex clustering [76j. 

MgB 2 is a ten year old new inter-metallic superconduc- 
tor which has many fascinating properties, particularly 
a high T c sa 39 K. We briefly discussed its electronic 
and geometrical structure in the pure case. For under- 
standing the occurrence of superconductivity we showed 
in Fig. 7 its two sets of a and 7r bands. The electrons 
in these bands are interacting via Coulomb and phonon 



interactions. There are unambiguous calculations and ex- 
periments which suggest the high strengths of inter-band 
and comparable intra-band interactions in this two-band 
material. It has been a consensus that MgB2 is a phonon- 
mediated superconductor. The higher T c is a result of 
coupling between the two Fermi surfaces (see e.g. Geerk 
et al. in |6Cf). This information provides a compelling 
reason to adopt a two-band model, either by microscopic 
BCS-Eliashberg or GL theory (see e.g. [52|, |53|, Zhito- 
mirsky and Dao in |55j, [77], and the review [78J])- Thus 
an argument of two independent bands is never tenable 
for MgB2 that exhibits one single T c . Hence no separate 
two superconducting regimes can occur in MgB2- 

A single vortex in a two-band superconductor will have 
a complicated vortex core that results from solving two 
coupled order parameter equations. Contrary to the re- 
mark by Babaev and Speight [60( , the core size and shape 
depends nontrivially on the coupling of two order param- 
eters and also on temperature. The next issue is about 
the penetration depth A. This quantity comes out of 
the GL equation for the vector potential. In the pres- 
ence of electromagnetic field coupling with the order pa- 
rameter [72 (n* - )/)^" 1X02) + he] as mentioned earlier in 
Sec. 4, A depends on the supcrfluid densities in a non- 
trivial fashion, see Eq. 9 and Fig. 3 of [79J . In brief, mix- 
ing two components as visualized by Babaev et al. and 
Moshchalkov et al., will not produce a new kind of super- 
conductivity, because in the mixture A is determined self- 
consistently. Table 1 of Golubov et al. [74| lists 12 val- 
ues of the zero-temperature London penetration depths 
Al(0) from various experiments with MgB2; these val- 
ues exceed the reported values of the coherence lengths 
£, and thus k ^> 0.71. Therefore, type-I condition will 
be hardly satisfied. For this reason it will be instructive 
to assess the "semi-Meissner" or "type- 1.5" conditions in 
relation to Moshchalkov's experiments [49j,|50|. 

Another question is: Can one have two different pene- 
tration depths A in the same material and for the same 
direction of the supercurrent density j ? Usually A is 
defined via the general linear relation (valid for small 
j and small B in isotropic superconductors, see Sec. 3) 
j(k) = -Q(fc)A(k) by Q(0) = A~ 2 or A = Q(O)- 1 / 2 . In 
general one has j = SF/SA where F{j, A} is the free- 
energy functional of a given theory and A the vector 
potential. This definition allows only for one single A. 

In conclusion, from our careful analysis we argue that 
the hypothesis of "type-1.5 superconductivity" as put for- 
ward in [49j,|50| and in the abstract of Moshchalkov's talk 
at this conference, is unrealistic and unfounded for MgB2 
that has one single transition temperature T c . Anyway, it 
cannot be precluded that the vortex-vortex interaction, 
if calculated microscopically (e.g., by the quasiclassical 
Eilenberger method [3l|, [32J, [8fJ, |81| ) may have an attrac- 
tive tail due to the complexity of the problem. 
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